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ABSTRACT 

We study the correlation heights, which indicate the formation height of Ex- 
treme Ultraviolet (EUV) lines in an active region using observations from the 
EUV Imaging Spectrometer (EIS) and Solar Optical Telescope (SOT) on board 
Hinode. The nonlinear force-free field (NLFFF) optimization method is adopted 
to extrapolate the 3D magnetic fields to higher layers. Three subregions with dif- 
ferent characteristics are selected in the active region for this study. The results 
show that the formation heights in different subregions vary with their different 
magnetic fields or velocity patterns. After solving the line blending problem be- 
tween the He II 256.32 A and Si x 256.37 A lines by the double Gaussian curve 
fitting, we find that the transition region lies higher in a strong magnetic area. 
In a pre-flare heating area there possibly exist multithermal loops as implied by 
comparing the Doppler velocity and the magnetic field on the solar disk. 

Subject headings: Sun: corona — Sun: magnetic fields — Sun: transition region 
— Sun: UV radiation 
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Introduction 



The correlation height, defined as the height at which the correlation coefficient between 
the Doppler velocity of extreme ultraviolet (EUV) lines and the magnetic field i nclination 



reache s its maximum, can be used to represent the formation altitude of EUV lines (ITu et al. 
2005al fbh. There are two reasons for this. On the one hand, the magnetic field controls the 



plasma motion in the transition region and corona, where the low (3 condition takes the 
responsibility. The plasma moves mainly along the magnetic field lines. Therefore the line- 
of-sight component of the velocity should be proportional to the cosine of the angle between 
the magnetic field lines and the line of sight, or the magnetic field inclination. Consequently, 
the absolute value of the Doppler velocity has a positive correlation with the value of \B Z /B\, 
where B is the magnetic field strength and B z is the line-of-sight component. Note that the 
cartesian coordinate system is adopted here. On the other hand, the EUV line intensity 
from an optically thin plasma is proportional to the integration j v G(T)n 2 e dV, where G(T) 
i s the contribution fun ction for plasma with temperature T and n e is the electron density 
(IWatanabe et al.l 120071 ). Although the volume V has a spatial extention, the contribution 
function G(T) peaks sharply at a temperature T max . The line emission is mostly from the 
plasma with temperature T max . Therefore, we can take the height where the temperature 
reaches T max as the EUV line formation altitude. 

The formation heights of EUV lines are very important for us to understand the struc- 
ture of th e solar upper atmosphere, wh ich has been simulated through 3D MH D numerical 



mode ling (IGudiksen fe Nordlundl 120051 ) and synthesis of coronal emission lines (IPeter et al. 
2004 ). However, it is d ifficult to deduce the formation height from observations. The first 



effort has been made by lTu et al.l ( 2005aT) who stud ied the transition region lines in the quiet 
Sun by the correl ation method. Tu et al.l (l2005bl ) further proposed the funnel model for 
solar wind origin. iMarsch et all tood ) studied the structure of the solar transition region 
in a polar coronal hole, and found that the correlation heights differ in the open magnetic 



field region and the closed on e even for the same emission line. IMarsch et al. 



(12004 l2008h 



and IWiegelmann et al.l (120051 ) studied the link between plasma flows and magnetic fields in 
active regions and in an equatorial coronal hole, respectively. The emission heights of coronal 
bright points derived from th e correlation of Fe XII intensity enhancement and t he magnetic 



fields have been discussed in iTian et al.l ( 120071 ). With the correlation method, iTian et al. 



(120081 ) studied the sizes of transition region structures in coroanl holes and the adjacent 
quiet Sun. Moreover, the magnetic field structure of qu iet Sun and coronal hol e s has been 
investigated with the help of a potential field model in IWiegelmann fe Solankil (12004 ) and 
compared with transition region and coronal EUV images. 



The correlation method mentioned above needs the combination of observed EUV lines 
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and 3D magnetic fields, the latter of which are usually obtained by force-free field extrap- 
olation of the measured fields in the lower solar atmosp here, such as the photosphere. In 

3hl l2002h 



or 



previous papers a linear force-free field (LFFF) model (I Wiegelmann fc Neukirc 
potential field model was adopted . Recently, a new nonlin ear force-free field ( NLFFF) opti- 
mizat ion method was proposed by lWheatland et al.l (120001 ) and implemented by IWiegelmann 
( 12004] ). which has been widely used for the extrapolation of 3D magnetic fields. The NLFFF 
extrapolatio n method is more accurate for magnetic field reconstruction especially in ac- 
tive regions (jSchrijver et al.ll2006l ; iMetcalf et al.l 120081 ) . where the magnetic fields are highly 
sheared and twisted, the electric current is large, and the force free parameter a usually 
changes through the space. The potential field and LFFF model are ina dequate for these 
active regions. A review on this topic can be found in I Wiegelmann! (120081 ). 



In this paper we study the correlation heights between the Doppler velocity of EUV lines 



Spectrometer (EIS: ( 


^ulhane et al. 


2007 


) and Solar Optical Telescoi 


2008 


Suematsu et al 


2008: 


Ichimoto et al. 


2008; 


Shimizu et al.ll2008l 



2007). In particular we adopt the NLFFF optimization method implemented by IWiegelmann 
( 120041 ) to extrapolate the 3D magnetic fields. Three subregions with different characteristics 
are selected in the active region for this study. The details of data analysis and image co- 
alignment are described in § [2j Results are presented in § [31 Discussion and conclusions are 
made in § HI 



2. Observations and Data Analysis 

2.1. EIS Doppler Velocity Fitting 

We select six EUV emission lines as shown in Table [T]to calculate the correlation heights. 
The six emission lines cover a relatively large range of formation temperatures, from the 
transition region to the corona. They were obtained by the EIS 1" x 256" slit scanning from 
west to east during 19:07 - 23:46 UT on 2006 December 12 in the active region NOAA 10930. 
The exposure time is 30 s for each slit. We apply the standard EIS procedure to calibrate the 
data and fit the lines with single Gaussian curve. The Doppler velocities can be calculated 
by the line center shifts. When determining the center of a specific line, two instrumental 
effects should be taken into account. One is the slit tilt, because the EIS slit is not perfectly 
parallel to the pixel array of the CCD. The other is the orbital variation, which is caused by 
thermal effect on the instrument and follows a sinusoidal behavior. To remove the former, 
we subtract the Doppler shift taken early in the mission (regarded as purely due to the tilt 
effect) from the current one. To remove the latter, we first set a rest line center position for 
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each raster which is taken to be the average of the southern part (a relatively quiet region) 
of each slit. Then the time variation of the rest line center relative to its theoretical position 
in laboratory is obtained. After the above steps, we can finally construct the Dopplergrams 
for each line. 

The average values of x 2 m the line profile fitting in some regions of interest, which will 
be defined in § 12.41 are shown in Table [TJ They are relatively small for Fe vill, Fe x, Fe 
xm, and Fe XV lines compared to He II and Fe XI lines in all the regions, since the latter 
two lines are blended with other lines. The He II 256.32 A line is blended with the Si x 
256.37 A line, especially in active regions and magnetic loops. Therefore, we use a double 
Gaussian curves to fit the two components with only one assumption that the blueward 
co mponent correspond s to the He II line. The Fe XI line is a doublet at 188.23 A and 188.30 



A (jYoung et al.ll2007al ). We also fit it by double Gaussian curves assuming that the 188.23 
A component lies blueward. Note that the double Gaussian curve fitting is sensitive to 
the selection of initial values, in particular the line center wavelengths. Here, we set them 
to be the rest wavelengths of each line. The fitting yields seven parameters, including the 
amplitudes, widths, and central wavelengths for both components, and a shared background. 
The result of double Gaussian curve fitting for the He II line is shown in Figure [T] and Table 
[TJ We find that the average velocities are more blue shifted compared to the result from 
single Gaussian curve fitting, as expected. In particular, the \ 2 values are much smaller 
in the former show in Table [T], implying that the fitting goodness is significantly 

improved. The result for the Fe XI 188.23 A line is similar to the He II line. 

The non-Gaus sian line profiles of EIS in this active region have been reported by 



Imada et al. ( 2008 ). who fitted the Fe xiv 274.20 A line in the arcades during the flare. 



Asai et al.l (j2008l ) performed the first spectroscopic observation of an MHD fast mode shock 
wave during the flare using double Gaussian curve fitting of the lines. Here, we analyze the 
lines before the flare. 



2.2. Magnetic Field Extrapolation 



The active region NOAA 10930 has been studied by many authors (IKubo et a" 



2007; 



Zhang et al.ll2007l : Ijing et al.ll2008l : ISchriiver et all2008l : Wang et all200a ) . iGuo et al.l (12008 ) 
extrapolated the 3D magnetic field using the NLFFF optimization method and studied the 
configuration, energy release and shear angle changes in this active region. The photospheric 
vector magnetic fields are obtained by fitting the polarized spectrum observed by SOT/SP 
(spectropolarimeter), which scanned the southwestern part of the solar disk from 20:30 to 
21:33 UT on 2006 December 12. In the fast map mode, the spatial sampling of SOT/SP 
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is about 0.3". The spectra are fitted by a nonlinear least-square fitting procedure based 
on the Milne-Eddington model. The 180° ambiguity of the azimuth angle is remoyed by 



comparing the obse rved field to the extrapolated LFFF field (jWangill997l ; IWang et al.l 12001 



Metcalf et al.ll2006l ). The top left panel in Figure [2] shows the full field of view (FOV) of 
the magnetic field observed by SOT/SP. We select half of the FOV and sample them to 
128 x 128 grids by 4 x 4 binning due to the limitation in comput ations. Then the boundar y 
data are preprocessed by a preprocessing routine developed by IWiegelmann et al.l (120 06). 
The final extrapol ation result by the NLFFF optimization method (IWheatland et al.ll2000l ; 
Wiegelmannl 120041 ) is shown in the bottom left panel of Figure [2j 



We adopt the potential field model to extrapolate the magnetic field in a larger FOV, 
because the extrapolation height of SOT/SP magnetic field is limited to ~ 110 Mm and 
the potential field model needs only the line-of-sight component of the field, which can be 
obtained in the MDI full disk observations. The FOV of MDI data used for extrapolation is 
shown in Figure [2] (top right), but the data in the area corresponding to the FOV of the top 
left panel are replaced by SOT/SP magnetic fields. This is because that the MD I data suffer 
from the saturation problem, which is discussed in detail in a previous paper (Moon et al. 



2007 . We extrapolate the fields with the Fourier transform method ( Alissandrakisl 



1981 



Gary! Il989l ) in this larger FOV that is four times the FOV used in the NLFFF method. 
The boundary data are sampled to 256 x 256 grids with spatial resolution of ~ 1.2". So 
the extrapolation height is ~ 220 Mm. The 3D magnetic field is shown in Figure [2] (bottom 
right). In order to calculate the correlation heights, we have to align the EUV images with 
the SOT/SP images. This is done by correlating the two set of images and finding the offsets 
between them. 



2.3. Co-alignment of Hinode Observations 



EIS covers two wavelength bands, i.e. short wavelength (SW) band (170 - 211 A) and 
long wavelength (LW) band (246 - 292 A), which are recorded by two CCDs, respectively. 
Images on the SW CCD have an offset of 16 - 20" northward and 2" westward compared 
to that on the LW CCD (lYoung et al.ll2007bl ). We correct for the offsets in SW images by 
applying a fixed value of 16" in the north-south direction and 2" in the east-west direction. 
We also need to align the SW and the LW images with the magnetograms. The main idea 
is to align all the images by EIS and SOT with that by EIT ( EUV Imaging Telesc ope; 
Delaboudiniere et al.l Il995l ) and MDI (Michelson Doppler Imager; IScherrer et al.l Il995l ) on 
SOHO (Solar and Heliospheric Observatory). The full disk images by EIT and MDI enable 
us to align them by comparing the limbs. The accuracy is close to their spatial resolutions. 
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Since the EIS and SOT observations are within partial disk, we cannot get the absolute 
coordinates of the solar disk center. However, after being aligned with the MDI images, all 
the coordinates are referred to the true solar disk center. 

The SOT/SP magnetogram is aligned with the MDI magnetogram as follows. First of 
all, the two images are interpolated to the same resolution, always to the better one. Then, 
we select two slit images to do the correlation because the MDI magnetogram was observed 
during the scanning time of the SOT/SP one, i.e. 20:30 - 21:33 UT. The slit image of 
SOT/SP that has the closest time to the MDI observation is correlated with the MDI slit 
images in different positions. Finally, the offset is found for which the two slit images have 
the largest correlation coefficient. The uncertainty of the co-alignment is about 2.0". 

The EIS Fe XII 195 A image is aligned with the EIT Fe XII 195 A image with the 
correlation method, too. Then the EIT 195 A image is aligned with the MDI magnetogram 
by comparing the limbs. Reading the coordinates of the east, west, south, and north poles on 
the solar limb in the MDI magnetogram, we calculate the coordinates of the disk center and 
find that the coordinates provided by MDI are correct. But the EIT 195 A image taken at 
23:49 UT on the same day, which is the closest to the EIS observation, shifts 4.0" westward 
and 6.0" northward in the coordinates provided by the standard software. The uncertainty 
of the shift is ~ 2.6", which is the resolution of EIT images. Finally, The offset between the 
EIS Fe XII 195 A image and the MDI image is found. The total uncertainty is ~ 5.2". 

The EIS He II 256 A velocity map is aligned with the SOT/BFI (Broadband Filter 
Imager) Ca II H image, the SOT/SP and MDI magnetograms consecutively. We find that 
the red shift region of the He II velocity map has similar configurations with the heating 
region of the Ca II H image as shown in Figure [TJ Note that the Ca II H image was observed 
every 2 min, but the EIS map was constructed from raster scanning during about 4.5 hours. 
We need to reconstruct a Ca II H image in the exactly same manner as the EIS map. To 
do so, we cut the Ca II H images at each time into slices along the slit. We pick out the 
slice with the same position and time as each slit in EIS scanning observations. Then we 
adopt a feature identification method to align them, i.e. select points with similar features 
alternately from the two images. Usually we select five points from each image to calculate 
an offset and repeat it 10 times to estimate the uncertainties. The average offsets in x and 
y directions are 4.3" and —7.8", with errors of 1.7" and 1.2", respectively. Next, the Ca 
II H image is aligned with the SOT/SP magnetogram with the identification method, too. 
The offset in x direction is —0.6" with an error of 0.4", and that in y direction is —7.1" 
with an error of 0.4". Because the SOP/SP magnetogram has already been aligned with the 
MDI magnetogram, the EIS He II velocity map can be aligned with the latter with total 
uncertainties of 4.1" and 3.6" along x and y directions, respectively. 
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We plot the offsets of all images with respec t to the MDI i mage i n Figure [3) If we adopt 
the offsets between S W and LW images given by lYoung et al.l (l2007bl ) , we only have to align 
one of them with MDI magnetograms. Then the other can be aligned automatically. But 
here we align both with independent methods to check the offsets derived by us and by them. 
The LW and SW images have the same offsets with respect to MDI magnetogram within 
the uncertainties, which indicates our results are coincident with theirs. 



2.4. Correlation Heights and Error Estimation 

We select different regions to calculate the correlation coefficients that are plotted in 
Figure HI Different from quiet Sun regions, the environment in active region is quite inho- 
mogeneous. It is necessary to select different subregions to calculate the correlation heights. 
The high spatial resolution observations of Hinode provide us the possibility to do so. The 
subregions are selected according to spatial patterns of the magnetic fields and velocity fields; 
so their shapes are not necessarily regular. Region 1 mainly concentrates on the negative 
magnetic polarity region and the blue shift area of He II 256 A. Region 2 is on the east side 
of the polarity inversion line corresponding to the foot points of some lower loops and the red 
shift area of He II 256 A. Region 3 is within the projection area of the higher closed magnetic 
loops. The correlation coefficient distributions are calculated by correlating the magnetic 
field inclination with the Doppler velocity, which is obtained by single (middle column) and 
double (right column) Gaussian curve fitting in Figure HJ respectively. The height is less 
than 5 Mm in Region 2, but more than 20 Mm in Regions 1 and 3. For different lines listed 
in Table HJ the correlation heights can be calculated with the same method. 

There are errors in the image co-alignment, which also cause errors in calculations of 
the correlation height. As mentioned in § I2.3[ there is an offset between SW and LW EIS 
images. Moreover, different lines in the same band have internal pointing shifts, which are 
0.5" in x direction and 1" in y direction for SW and 3" in x direction and 4" in y direction 
for LW as shown in a preliminary research by Det0. Taking into account these facts, the 
errors in the image co-alignment are 7.7" in x direction and 8.2" in y direction for SW except 
for Fe XII 195 A. They are 7.1" in x direction and 7.6" in y direction for LW except for He 
II 256 A. The errors of the correlation height are estimated by the following way. First, we 
artificially shift the EIS velocity map relative to the magnetic field B z image within the range 
of co-alignment errors. Then, the correlation heights are calculated for each case (here, we 
select 1" x 1" grids). Finally, we get a grid of correlation heights with different correlation 



1 http: / /msslxr.mssl.ucl.ac.uk:8080/eiswiki/Wiki.jsp?page=FitsPointing 
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coefficients. The final correlation height for one specific line is estimated as the average 
of the heights with correlation coefficients greater than 0.07, and the error is the standard 
deviation of these heights. 

3. Results 

The average formation height of He II 256 A in region 2 is lower than that in Regions 1 
and 3. In comparison, the mean magnetic fields in Region 2 and 3 are one order of magnitude 
weaker than that in Region 1. This implies that stronger magnetic fields correspond to higher 
formation heights of EUV lines. For Regions 2 and 3, the mean magnetic fields are almost 
the same, but the Doppler velocity patterns are different, which show mainly red shifts and 
blue shifts, respectively. This means that down flows tend to make the formation height 
lower. In summary, the correlation heights are different in various regions even for the same 
line He II 256 A because of different physical conditions. 

The correlation height versus log T max is plotted in Figure where the EUV line forma- 
tion temperature, T max , is listed in Table [1] and the errors of formation heights are derived 
using the method mentioned in § 12.41 In Region 1, shown in Figure [5] (top left), there are no 
points for Fe vm and Fe XV, since we cannot find the heights with correlation coefficients 
greater than 0.07. A possible reason is that the real formation height of Fe XV exceeds the 
NLFFF extrapolation domain, which is about 110 Mm in height. For Fe vm we still do not 
know the reason. The correlation heights for all the EUV lines can be got in Region 1 by 
the potential field extrapolation in Figure [5] (top right), which can reach about 220 Mm in 
height. For the regions studied here, the heights for the same line calculated with different 
extrapolation methods differ slightly within the errors. Also, the heights calculated by the 
potential field give a similar height versus log T max relationship compared to that calculated 
by NLFFF. However, the NLFFF gives less uncertainties for the correlation heights and 
makes the results more confirmative. 

Region 2, shown in Figure [5] (bottom left), is a pre-flare heating area. Massive downflows 
can be seen in the EIS velocity map (Figured]). All the EUV lines are formed within a smaller 
range of heights. The temperature structure is complicated in this area. There is no clear 
height versus log T max relationship in Region 2, while there are positive relationships in 
Region 1 and Region 3 (bottom right) as shown in Figure [5]. But the temperature in Region 
3 rises more quickly. For example, log T max = 6.3 corresponds a height of ~ 90 Mm in 
Region 3, compared to > 100 Mm in Region 1. This is probably due to the higher magnetic 
field strength in Region 1, which seems to suppress the temperature increase with height. 
Therefore, the EUV line formation height is affected by not only the magnetic topology but 



- 9- 



also the field strength. 

We also plot the linear fitting of the height - log T max relationship for Regions 1 and 3 
in Figure [5J These curves provide a quantitative description of how the temperature varies 
with height in the corona. 



Discussion and Conclusions 



The EUV lines in an active region are formed at different heights according to the 
physical conditions of the magnetized atmosphere in situ, including the velocity, density, 
temperature, and magnetic field. In this paper we examine some of these factors with the 
correlation height method by aligning images observed by SOT and EIS. In particular, we 
find that the transition region, where He II 256 A is formed, lies higher in areas with strong 
magnetic fields. 

Results from Region 2 show that the EUV emission lines are formed within a relatively 
narrow height range in the pre-flare heating l oops. A possible scenario is tha t this region 
contains coronal loops that are multithermal. iKjeldseth-Moe fc Brekkd (119981 ) have shown 



that some active region loops contain mixtures of plasma of differ ent tern 



peratures. Multi- 



thermal plasma has also been found in the transequatorial loops (iBrosiusi 120061 ) . All these 



researches are based on the analysis of the co-spatial emission lines on the solar limb. While 
we find a similar result by analyzing the lines on the solar disk and correlating the Doppler 
velocity and magnetic field. 
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Table 1: EIS emission lines, average Doppler velocities, and average values of x 2 for the line 
profile fitting in three selected regions. 



Ion Wavelength log T max Average Doppler Velocity a Average x 2a 

(A) (km s- 1 ) 

He ll b 256.32 4~9 ^8 43 ^3 L6 jfo 3^3 

He ll c -12.9 1.7 -2.7 0.5 0.4 0.5 

Fevm 185.21 5.6 -21.9 1.6 -5.2 0.3 0.3 0.3 

Fe x 184.54 6.0 -27.5 1.2 -13.8 0.4 0.7 0.6 

Fe xi b 188.23 6.1 -15.4 -1.0 -9.9 2.5 12.4 10.1 

Fe xi c -25.3 -9.7 -14.4 0.5 0.3 0.5 

Fe xm 202.04 6.2 -24.3 -4.0 -10.9 0.3 0.4 0.4 

Fe XV 284.16 6.3 -11.8 -5.3 -1.7 0.7 2.0 1.2 



a The values in the three columns are calculated for Regions 1-3, which are defined in Figure HI 
b The average velocities and \ 2 f° r this line are calculated by single Gaussian curve fitting. 
c The average velocities and x 2 for this line are calculated by double Gaussian curve fitting. 
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Fig. 1. — EIS He II 256 A velocity maps fitted by single Gaussian curve (left panel) and 
double Gaussian curves (middle panel). The velocities that excess 40/— 40 km s -1 have been 
set to the top color (Red/Blue) in order to make the velocity field pattern stand out. The 
Ca II H image (right panel) shows that it has similar configurations with the He II 256 A 
velocity maps. The Ca II H image observed by SOT/BFI is created artificially to match the 
EIS map in time. See text for details. 
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Fig. 2. — Top left: The FOV of the magnetic field observed by SOT/SP from 20:30 to 
21:33 UT on 2006 December 12. The dashed square denotes the FOV in which the NLFFF 
extrapolation is performed. Bottom left: The 3D magnetic field extrapolated by the NLFFF 
model for the dashed square shown in the top left panel. Top right: The FOV of MDI 
magnetogram for potential field extrapolation observed at 20:51 UT on 2006 December 12. 
The dashed square is the same to that in top left panel. Bottom right: The 3D magnetic field 
extrapolated by the potential field model. Note that only the region in the dashed square is 
shown. 
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Fig. 3. — Offsets of images observed by different instruments with respect to MDI. Plus 
values in the parentheses refer to that image shifts to east or south, while negative values 
to west or north. The error bars are also plotted, both in the east-west direction and the 
north-south direction. 
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Fig. 4. — Height distributions of the correlation coefficient for the He II 256 A line in 
different subregions of the active region NOAA 10930. The white polygons labelled with 
Regions 1-3 (left column) show the areas selected to calculate the correlation coefficients. 
The distributions in middle and right columns are calculated by correlating the magnetic 
field inclination with the Doppler velocity, which is obtained by single and double Gaussian 
curve fitting, respectively. The distributions in top to bottom rows are obtained in Regions 
1-3, respectively. 



150 



100 



50 



Region 1 : NLFFF 

h=a+b log(7' mo< ) 
a=-148.4±151.9 
6= 39. 6± 25.2 



[3-- 





4.5 

1 50 



5.0 



5.5 



100 



50 



Region 2 : NLFFF 




4.5 



5.0 



5.5 

Log J n 



6.0 



150 



1 00 



50 







Region 1 : Potential Field 
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a=-285.8±1 18.6 
6= 63. 3± 23.0 
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Region 3 : NLFFF 

h=a+b \oq(T m ) 
a=-223.5± 54.3 
6= 49. 8± 10.5 
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Fig. 5. — Correlation heights versus log T max for Regions 1-3. The magnetic field is 
extrapolated with the NLFFF model except for the top right panel, in which the potential 
field model is used. The dotted lines are the linear fitting of the height - log T max relationship. 



